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Rate  constants  for  the  reactions  of  02~  and  O'  with  N  and  O  atoms  have  been  measured  for  the 
first  time  as  a  function  of  temperature  from  173  to  500  K  for  O'  reactions  and  173  to  400  K  for 
O2'  reactions.  Room  temperature  rate  constants  for  (V  reacting  with  N  and  O  are  3.1  x  1 0  10  and 
1.7  x  10  10  cm3  s'1,  respectively,  and  the  corresponding  O  rate  constants  are  1.7  x  1 0  1,1  and 
1.5  x  10  10  cm3  s  1 .  in  good  agreement  with  previous  values.  Temperature  dependences  are  about 
T'1,7  for  both  CT  -  reactions  and  T  06  and  T  1,3  for  the  reactions  of  O'  with  N  and  O,  respec¬ 
tively.  Branching  for  the  CB”  reaction  with  N  is  found  to  predominantly  form  O-  (>85%)  in 
contrast  to  previous  measurements,  which  reported  NCB  +  e~  as  the  main  channel.  Calculations 
point  to  the  present  results  being  correct.  The  potential  energy  surface  for  this  reaction  was  calcu¬ 
lated  using  density  functional  theory,  coupled  cluster  with  singles,  doubles  (triples),  complete  ac¬ 
tive  space  self-consistent  held,  and  complete  active  space  second-order  perturbation  methods  and 
is  found  to  be  quite  complex,  with  agreement  between  the  calculated  surface  and  the  observed  ki¬ 
netic  data  only  possible  through  the  inclusion  of  dynamical  correlation.  ©  2013  AIP  Publishing  LLC. 
[http://dx.doi.org/ 10.1 063/1 .48240 1 8] 


INTRODUCTION 

Already  in  1961,  Dalgarno  suggested  that  O  and  CB” 
react  with  O  atoms  to  detach  electrons.1  This  was  later  con¬ 
firmed  by  the  NOAA  group  in  1967,  who  also  found  that  N 
atoms  detach  electrons  in  reaction  with  O'  and  O2-.2  The 
potential  pathways  are  listed  in  reactions  (l)-(4).  The  NOAA 
group  reported  that  reaction  (la)  dominated.  They  further 
established  that  (2a)  and  (2b)  competed  but  did  not  report 
branching  ratios.  These  findings  proved  especially  important 
for  the  chemistry  of  the  D-region  of  the  ionosphere  because 
the  reactions  convert  heavy  ions  into  light  e~  which  impact 
radiowave  propagation:3  5 


r  +  N-* 

N02  +  e~, 

ArH°  = 

-4.11  eV, 

(la) 

NO  +  O', 

ArH°  = 

-2.39  eV, 

(lb) 

h  +  O  ~ 

►  03  +  e~, 

ArH°  = 

-0.66  eV, 

(2a) 

O2  +  O  , 

ArH°  = 

-1.01  eV, 

(2b) 

O'  +N 

— >  NO  +  e~ 

,  ArH° 

=  -5.09  eV, 

(3) 

0"-t-0 

-*  02  +  e“ 

,  ArH° 

=  -3.71  eV. 

(4) 

“^Author  to  whom  correspondence  should  be  addressed.  Electronic  mail: 
afrl.rvborgmailbox@kirtland.af.mil 


No  subsequent  work  had  been  performed  on  reactions 
( 1 )— (4)  until  our  group  reexamined  reactions  (1)  and  (2)  in 
2006, 6  as  part  of  an  ongoing  effort  to  improve  on  the  knowl¬ 
edge  of  ionospheric  chemistry.7  In  that  study,  we  found  that 
chemical  reaction  (lb)  competes  with  detachment  in  the  reac¬ 
tion  of  Oo-  with  N.  This  discrepancy  with  the  literature  was 
attributed  to  using  a  selected  ion  flow  tube  with  an  external 
ion  source  rather  than  a  flowing  afterglow.  Here,  we  refine 
the  measurements  and  report  the  first  temperature  dependent 
measurements  for  all  four  reactions,  and  detail  our  modeling 
efforts  for  reaction  (1). 


EXPERIMENTAL 

All  experiments  were  performed  on  a  selected-ion  flow 
tube  instrument  that  has  been  described  in  detail  in  Refs.  6,  8, 
and  9.  Briefly,  CB-  was  generated  by  associative  electron  at¬ 
tachment  to  O2,  while  O'  was  created  by  dissociative  attach¬ 
ment  to  N2O  in  a  high-pressure  electron  impact  ionization 
source.  The  ions  were  focused  into  a  quadrupole  mass  filter, 
mass  selected,  and  injected  into  the  flow  reactor  via  a  Venturi- 
type  inlet.  Ions  undergo  ~105  collisions  with  helium  buffer 
gas  at  pressures  of  0.4  to  0.8  Torr  resulting  in  complete  or 
near-complete  thermalization. 111  The  higher  pressure  was  used 
when  studying  the  high  temperature  reactions  because  the  rate 
constants  are  small  and  only  limited  amounts  of  atoms  were 
formed.  Higher  pressure  was  achieved  by  throttling  a  butterfly 
gate  valve  resulting  in  lower  pumping  speeds  and  thus  longer 
reaction  times.  Neutrals  were  injected  49  cm  before  the  end 
of  the  flow  tube  and  allowed  to  react  with  02“  or  O'.  The 
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product  ions  and  residual  reactant  ions  were  sampled  through 
a  1  mm  orifice  into  a  second  quadrupole  mass  filter  for  mass 
analysis.  A  small  attractive  voltage  (~5  V)  was  applied  to  the 
nose  cone.  Ions  were  detected  with  an  electron  multiplier  op¬ 
erating  in  pulse  counting  mode.  The  temperature  of  the  flow 
tube  was  varied  from  173  to  500  K  either  by  resistive  heating 
devices  or  pulsed  liquid  nitrogen. 

Nitrogen  atoms  were  generated  via  microwave  discharge 
of  N2  in  a  Pyrex  side  arm.  The  extent  of  nitrogen  dissoci¬ 
ation  was  small  (<1%)  and  was  varied  by  changing  either 
the  N2  flow  rate  or  the  microwave  discharge  power.  Contri¬ 
butions  from  vibrationally  excited  N2  from  this  setup  have 
been  shown  to  be  negligible  due  to  collisions  in  the  glass 
side  arm. 11,12  Oxygen  atoms  were  generated  using  a  titration 
method  that  has  been  described  in  detail  in  Refs.  11  and  13. 
The  surface  of  the  glass  was  treated  with  a  boric  acid  solution 
to  minimize  O  atom  losses  on  the  walls  of  the  cavity.  NO  gas 
was  introduced  at  a  known  flow  rate  into  the  glass  side  arm 
downstream  of  the  discharge  where  it  quantitatively  reacted 
with  the  nitrogen  atoms  to  give  ground  state  oxygen  atoms.11 
A  run  started  by  recording  ion  signals  with  the  discharge  off 
then  on  without  addition  of  NO.  The  decrease  in  the  signal 
intensity  combined  with  the  N  atom  density  determined  the  N 
atom  rate  constant.  NO  was  then  added,  which  converted  the 
N  into  O.  At  low  NO  flows,  both  N  and  O  atoms  flowed  into 
the  reaction  zone  and  the  NO  was  completely  converted  to  O. 
At  higher  flows,  O  and  NO  both  entered  the  reaction  zone.  A 
change  of  slope  occurred  at  the  point  where  all  N  atoms  were 
converted  to  O  atoms,  i.e.,  the  titration  end  point. 


RESULTS  AND  DISCUSSION 

Figure  1  shows  a  set  of  raw  data  for  300  K  and  0.8 
Torr  for  the  02~  reactions.  The  black  triangle  represents  the 
02~  signal  with  the  discharge  off;  the  O-  signal  without 
the  discharge  is  near  zero  and  is  not  displayed.  Solid  circles 
and  squares  represent  discharge-on  data  for  O2  and  O  ,  re¬ 
spectively.  Turning  the  discharge  on  creates  N  atoms  which 


FIG.  1.  02~  and  O-  ion  signals  as  a  function  of  NO  flow  at  300  K  and 
0.8  Torr  after  a  reaction  time  of  5.8  ms.  The  black  triangle  denotes  data  taken 
before  the  discharge  was  turned  on. 


decrease  the  02_  signal  substantially  while  increasing  the  O 
signal  due  to  reaction  (lb).  As  NO  was  added  in  low  con¬ 
centrations,  converting  some  of  the  N  to  O,  the  ()2~  signal 
increased  indicating  the  O  atom  rate  constant  is  smaller  than 
the  N  atom  rate  constant  with  02  .  At  NO  concentrations  ex¬ 
ceeding  the  initial  N  atom  concentration,  the  amount  of  O  re¬ 
mained  constant  with  increasing  amounts  of  NO  added.  Since 
02_  does  not  react  with  NO,  the  signal  remains  constant  as 
expected.  The  NO  concentration  at  which  the  slope  changes 
from  positive  to  zero  equals  the  initial  N  atom  concentration, 
[N]().  The  O"  similarly  increased  at  low  NO  addition  but  then, 
because  O  reacts  with  NO,  decreased  after  the  titration  end 
point.  The  rate  constant  of  02_  with  N  is  derived  from  the 
following  equation: 


in  IQ2  W/ 

[O2  ]o ,on 

=  [N  ]0t 


(5) 


where  the  02_  signals  are  with  the  discharge  off  and  on  at 
zero  NO  flow  and  t  is  the  reaction  time.  The  positive  slope  at 
low  NO  concentrations  reflects  the  amount  that  the  O  atom 
rate  constant  is  smaller  than  the  N  atom  rate  constant.  The  O 
atom  rate  constant  can  be  derived  similarly  to  (5),  replacing 
the  discharge-on  02_  signal  at  zero  NO  flow  with  the  signal  at 
the  titration  end  point,  where  the  O  atom  concentration  equals 
[N]o-  Separate  experiments  were  performed  for  determining 
the  O-  rate  constants,  in  which  O-  was  injected  into  the  flow 
tube  instead  of  02_.  Absolute  uncertainties  in  the  rate  con¬ 
stants  are  on  the  order  of  ±40%  and  relative  rate  constants 
have  uncertainties  of  ±25%.  These  are  larger  than  the  uncer¬ 
tainties  normally  reported  for  SIFT  measurements  of  ±25% 
and  ±15%,  respectively.  The  larger  error  comes  from  more 
scatter  in  the  data  and  the  potential  for  atom  loss  on  the  flow 
tube  walls.  The  latter  effect  cannot  be  measured  so  we  treat 
these  errors  as  estimates  only.  Relative  errors  (from  one  tem¬ 
perature  to  the  next)  are  smaller  since  many  of  the  errors  are 
the  same.  Each  reported  rate  constant  is  the  average  of  about 
four  runs. 

Table  I  lists  the  present  results  along  with  literature  val¬ 
ues  of  the  rate  constants  for  300  K.  In  general,  all  values  are 
in  good  agreement.  The  present  results  are  about  25%  lower 
than  those  of  Fehsenfeld  et  al.2  with  the  exception  of  the  CF- 
with  O  reaction.  In  this  case  the  two  values  disagree  by  a  fac¬ 
tor  of  1.9,  which  still  is  within  the  combined  experimental 
error.  Our  agreement  with  the  previous  value  from  our  labo¬ 
ratory  (Poutsma  et  al.6)  for  the  02_  with  N  reaction  is  sim¬ 
ilarly  good,  although  in  the  opposite  direction.  However,  the 


TABLE  I.  Comparison  of  the  reaction  rate  constants  for  reactions  (l)-(4) 
measured  at  300  K  in  units  of  10-10  cm3  s-1.  The  last  column  is  the  temper¬ 
ature  dependence  measured  in  the  current  study  (expressed  as  fits  to  a  power 
law). 


Reaction 

Present 

Poutsma  et  al.6 

Fehsenfeld  et  al.2 

n  in  T-n 

02~  +N 

3.1 

2.3 

4.0 

1.7  ±  0.3 

O2  +  0 

1.7 

3.9 

3.3 

1.8  ±  0.3 

O-  +N 

1.7 

2.2 

0.6  ±  0.3 

O-  +0 

1.5 

1.9 

1.3  ±  0.3 
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present  CL”  with  O  rate  constant  is  over  a  factor  of  two  lower 
than  our  previously  reported  value.  While  the  two  previous 
values  for  O2  with  O  (i.e.,  Poutsma  et  al.  and  Fehsenfeld 
et  al.)  appear  to  be  in  good  agreement  with  each  other,  they 
have  an  opposite  ordering  of  the  N  and  O  reactions,  indicating 
the  difficulty  of  the  measurements.  We  speculate  that  our  cur¬ 
rent  value  is  more  accurate  for  several  reasons.  Looking  at  the 
raw  data  figures  (see  Fig.  1  in  both  papers)  shows  a  leveling 
out  at  high  NO  flows  in  the  current  study  as  expected,  since 
NO  does  not  react  with  CL',  while  the  previous  study  shows 
a  slight  decrease  in  the  CL'  signal  after  the  titration  end  point. 
This  indicates  that  the  previous  experiment  had  some  signal 
drift  as  O  atoms  were  added,  presumably  due  to  surface  ox¬ 
idation  causing  non-conducting  layers.  In  the  past,  we  were 
not  able  to  measure  temperature  dependences  for  this  reac¬ 
tion,  because  signal  drift  was  worse  than  in  our  current  setup, 
particularly  above  and  below  room  temperature.  A  key  dif¬ 
ference  between  the  instrument  as  employed  in  our  previous 
study  (at  Hanscom  AFB)  and  that  at  present  is  that  all  vac¬ 
uum  pumps  have  been  replaced  with  oil-free  pumps;  a  cleaner 
system  may  explain  the  difference. 

Figure  2  shows  the  temperature  dependences  of  the  four 
reactions.  Reactions  were  studied  from  173  K  to  400  K  for 
CL'  and  173  K  to  500  K  for  O  .  The  reduced  upper  tempera¬ 
ture  for  the  CL“  reactions  was  due  to  the  presence  of  thermal 
detachment  for  O2  at  500  K:3 


O-,  -f-  Fie  — >  CL 


-r  ..v  02  +  e"  +  He  ArH°  =  0.448  eV.  (6) 


Evidence  for  reaction  (6)  is  presented  below.  All  data  show 
a  clear  trend  that  increasing  temperature  decreases  rate 
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constants,  as  is  often  the  case  for  ion-molecule  reactions.14' 15 
The  fits  shown  in  Fig.  2  are  power  laws  and  the  exponents 
are  listed  in  Table  I.  Both  CL”  reactions  have  similar  temper¬ 
ature  dependences,  T”1,7  and  ’I’  l  x  for  the  N  and  O  reactions, 
respectively.  The  O  reactions  have  smaller  temperature  de¬ 
pendences,  T'0,6  and  T”1 3  for  the  N  and  O  reactions,  respec¬ 
tively.  The  smaller  values  for  O  reactions  are  expected  since 
the  reactants  have  no  internal  degrees  of  freedom. 

Branching  ratios  for  the  reaction  of  O2”  +  N  were  ob¬ 
tained  by  monitoring  the  intensity  of  both  the  O2”  reactant 
and  O  product  as  a  function  of  extent  of  reaction.  The  N 
concentration  was  varied  by  changing  either  the  microwave 
discharge  power  or  the  N2  flow.  Figure  3  shows  the  raw  data 
where  the  y-axis  is  the  O'  count  increase  divided  by  the 
CL”  count  decrease  (the  CL”  decrease  is  equal  to  the  sum 
of  all  products  produced).  The  x-axis  is  the  extent  of  reac¬ 


tion,  In 


[O2  ]off 

[02"U 


,  where  [O2  ]0ff  is  the  signal  with  the  discharge 


off.  There  is  a  decline  in  the  O'  fraction  with  increasing  ex¬ 
tent  of  reaction  due  to  reaction  (3).  A  linear  least  squares  fit 
indicates  100%  (from  the  y-intercept  at  110%)  branching  to 
produce  O'  (reaction  (lb)).  Looking  at  the  graph  we  can  give 
a  lower  limit  of  85%  O'  production.  Previously,  we  have  re¬ 
ported  this  as  only  35%. 6  The  difference  should  only  be  due 
to  mass  discrimination.  In  between  the  two  experiments  we 
upgraded  our  quadrupole  control  electronics  so  that  a  much 
higher  frequency  supply  is  now  available  for  experiments  in 
this  mass  range,  2.1  MHz  vs.  1.2  MHz.  We  measured  mass 
discrimination  in  the  current  study  by  studying  the  reaction  of 
0+  with  O2  producing  O24",  i.e.,  the  same  masses.  Discrim¬ 
ination  was  found  to  be  less  than  6%.  Obviously,  the  same 
test  cannot  be  done  for  the  old  measurements,  but  discrimina¬ 
tion  at  low  masses  was  sometimes  a  problem  with  the  lower 
frequency  electronics,  though  the  extent  of  it  appears  to  have 
been  missed  in  this  case  because  the  O'  product  also  detaches 
electrons.  At  temperatures  of  400  K  and  below  we  found  that 
within  our  uncertainty  (~20%),  only  O'  is  produced,  i.e., 
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FIG.  2.  Rate  constants  for  the  reactions  of  O2  and  O  with  N  and  O  as  a 
function  of  temperature. 


FIG.  3.  Data  used  to  derive  branching  for  reaction  (1).  The  y-axis  is  the 
amount  of  O  increase  divided  by  the  O2'  decrease  ( )  and  the  x-axis  is 
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reaction  (lb)  dominates  reaction  (lb).  At  500  K,  we  measured 
more  O-  produced  than  Oi  declined.  We  attribute  that  to  the 
presence  of  reaction  (6),  i.e.,  the  N  atoms  reacted  with  some 
of  the  O2-  before  it  detached.  Modeling  the  kinetics  con¬ 
firms  this  possibility  but  deriving  accurate  rate  constants  for 
reaction  (6)  is  not  possible.  The  branching  ratios  for  reaction 
of  C>2~  +  O  were  not  determined  in  the  present  work  due  to 
both  the  presence  of  either  N  atoms  or  NO  along  with  the  O 
atoms,  as  well  as  secondary  chemistry  greatly  complicating 
the  analysis  and  preventing  any  meaningful  results. 

These  measurements  have  direct  relevance  to  meso¬ 
spheric  chemistry.  Both  O  and  CF-  are  found  in  apprecia¬ 
ble  abundance,  as  well  as  O  atoms  and  to  a  lesser  extent  N 
atoms.3, 16-18  Since  electrons  interact  with  radio  waves  much 
more  extensively  than  ions,  processes  that  convert  ions  to 
electrons  are  important  for  communications.  Our  finding  that 
N  atoms  do  not  directly  produce  electrons  in  reaction  (1), 
but  instead  indirectly  since  the  O  produced  reacts  by  de¬ 
taching  an  electron,  will  not  have  a  major  effect  since  the 
N  atom  concentration  in  the  mesosphere  is  much  lower  than 
that  of  O.  What  will  have  a  larger  effect  is  that  the  rate  con¬ 
stants  are  all  larger  at  low  temperature  and  the  mesosphere 
is  typically  cold,  which  should  shift  ions  to  electrons  more 
rapidly.  Insight  into  reactions  (2) — (4)  is  difficult,  as  these  re¬ 
actions  are  challenging  to  model.  The  atom  reactions  involve 
electron  reactions,  cannot  be  modeled  statistically,  and  are 
highly  sensitive  to  details  of  the  potential  energy  surfaces.  The 
C>2~  reaction  with  O  is  charge  transfer  producing  O2,  which 
is  notoriously  difficult  to  model  properly.11’  More  tractable 
theoretically  is  reaction  (1),  for  which  we  have  calculated 
the  potential  energy  surface.  The  electronic  state  responsi¬ 
ble  for  the  reaction  was  found  to  be  the  lowest  triplet  (13A") 
state.  There  have  been  very  few  ab  initio  studies  of  the  NO2 
anion  on  this  triplet  state  and  the  existing  studies20,21  only 
focused  on  the  equilibrium  geometry  and  energy  without  ad¬ 
dressing  the  N(4S)  +  02-(X2ng)  reaction.  As  a  result,  the 
reaction  pathway  leading  to  the  formation  of  NCF-  is  sur¬ 
veyed  here  using  several  different  theoretical  methods.  Two 
minima,  both  with  C2V  symmetry,  have  been  located  at  several 
different  levels  of  theory,  and  their  geometries,  energies,  and 
vibrational  frequencies  are  summarized  in  Fig.  4.  The  global 
minimum  features  a  larger  O-N-O  bond  angle  than  the  higher 
local  minimum,  and  there  exists  an  isomerization  barrier  be¬ 
tween  the  two,  which  is  not  determined  here.  The  situation 
here  is  quite  similar  to  03(1 3 A"), 22  which  is  isoelectronic  to 
N02_(13A").  Density  functional  theory  (DFT)  calculations 
with  two  different  functionals  (MPW1K23  and  ®B97X-D24) 
give  similar  results,  while  the  spin-unrestricted  explicitly  cor¬ 
related  coupled  cluster  with  singles,  doubles,  and  perturba¬ 
tive  triples  calculations  (UCCSD(T)-F12)2:’  reduce  the  depth 
of  the  two  minima  slightly.  These  calculations  used  the  aug¬ 
mented  correlation-consistent  polarized  triple-zeta  (AVTZ) 
basis  set.26  The  UCCSD(T)  and  DFT  calculations  were 
performed  with  MOLPRO27  and  GAUSSIAN  09. 28 

Unfortunately,  single  reference  approaches  do  not  con¬ 
verge  well  in  the  reactant  asymptote  and  other  regions.  The 
three-dimensional  potential  energy  surface  mapped  at  the 
DFT  level  is  quite  bumpy,  making  it  difficult  to  identify  the  re¬ 
action  pathway.  The  difficulties  associated  with  the  electronic 
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FIG.  4.  Stationary  points  on  the  NCF  anion  (13A"  state)  potential  energy 
surface,  determined  by  various  theoretical  methods,  as  indicated  in  the  upper 
right  comer.  The  two  minima  are  both  in  C2v  symmetry,  and  the  correspond¬ 
ing  bond  lengths  and  bond  angles  are  given  in  angstroms  and  degrees,  respec¬ 
tively.  Energies  are  listed  in  kcal/mol,  while  frequencies  are  given  in  cm~ 1 . 
Blue  and  red  balls  represent  N  and  O  atoms,  respectively.  The  vibrational 
frequencies  are  listed  in  increasing  order. 


structure  calculations  are  presumably  due  to  the  fact  that 
there  are  many  low-lying  electronic  states  in  the  dissociation 
asymptote,  even  near  the  equilibrium,  as  was  found  for  the  O3 
molecule.22  As  a  result,  a  multi-reference  method,  namely, 
the  complete  active  space  self-consistent  field  (CASSCF) 
method20  implemented  in  MOLPRO27  was  used.  It  was  noted 
that  the  valence  active  space  typically  failed  to  converge,  but 
the  inclusion  of  the  3s  orbitals  of  the  O  atoms  was  found 
to  improve  convergence.  In  order  to  reduce  the  computa¬ 
tional  costs,  the  Is  and  2s  orbitals  of  N  and  O  were  all 
treated  as  doubly  occupied  but  the  molecular  orbital  coef¬ 
ficients  were  fully  optimized.  The  resulting  active  space  of 
(1  lo,  12e)  yielded  49104  configuration  state  functions.  In  ad¬ 
dition,  state-averaged  CASSCF  wave  functions  which  include 
the  two  lowest  A"  states  with  equal  weights  were  employed. 

As  shown  in  Fig.  4,  the  optimized  minima  at  the  CASSCF 
level  have  similar  geometries  as  those  obtained  from  other 
methods,  but  have  much  smaller  well  depths  relative  to  the 
CF-  +  N  reactant  asymptote.  This  is  apparently  due  to  the  ne¬ 
glect  of  dynamical  correlation  energy  in  the  CASSCF  method. 
Convergence  remains  an  issue  in  the  CASSCF  calculations 
in  the  dissociation  asymptote.  After  many  attempts,  we  ob¬ 
tained  a  reaction  path  along  the  Jacobi  coordinate  R ,  which 
is  the  distance  between  N  and  the  center  of  mass  of  C>2_, 
with  other  degrees  of  freedom  relaxed,  employing  the  rela¬ 
tively  small  VDZ  basis  set.  As  shown  in  Fig.  5,  it  has  a  small 
barrier  (dNo  =  1.9  A,  doo'  =  1-36  A,  0qno'  —  116°)  about 
0.6  kcal/mol  above  the  reactants.  It  is  interesting,  however,  the 
barrier  disappears  at  the  complete  active  space  second-order 
perturbation  (denoted  as  CASPT2  in  MOLPRO)  level,30  which 
captures  part  of  dynamical  correlation.  It  should  be  noted  that 
the  total  number  of  contracted  configuration  state  functions 
reached  84  687  237  in  the  CASPT2,  making  the  computa¬ 
tions  extremely  time-consuming.  Consequently,  higher  levels 
of  ab  initio  methods,  such  as  the  multi-reference  configuration 
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FIG.  5.  Potential  energy  curves  as  a  function  of  N-02~  Jacobi  distance  with 
other  coordinates  optimized  at  the  CASSCF/VDZ  level.  The  open  and  solid 
circles  correspond  to  the  CASSCF  and  CASPT2  energies,  respectively,  at  the 
CASSCF  optimized  geometries. 

interaction  (MRCI)  calculations,  are  not  presently  feasible  for 
this  molecule  due  to  formidable  computational  costs. 

The  negative  temperature  dependence  of  the  rate  constant 
suggests  a  barrierless  reaction  pathway.  Indeed,  the  CASPT2 
curve  presents  no  barrier,  consistent  with  the  experimental 
data  and  demonstrates  the  importance  of  dynamical  correla¬ 
tion  energy  in  this  system.  The  situation  here  is  similar  to  the 
case  of  O  +  O2  where  a  similar  entrance  channel  barrier  has 
been  shown  to  disappear  when  progressively  larger  basis  sets 
and  higher-level  ab  initio  methods  are  used.3 1  A  more  reliable 
determination  of  the  entrance  channel  potential  energy  sur¬ 
face  would  require  a  larger  basis  set  and  an  even  higher  level 
of  ab  initio  treatment.  Such  calculations  are  still  extremely 
challenging. 

The  observation  of  CV-  +  N  producing  primarily  NO  and 
O  ,  approaching  100%,  as  opposed  to  NO2  +  e~  deserves 
comment,  as  the  latter  channel  is  1.72  eV  more  exothermic 
than  the  former.  This  is  potentially  due  to  the  behavior  of 
the  specific  rate  curves  for  detachment  and  dissociation  from 
N02~.  As  seen  in  the  case  of  SFg-32  (typical  of  barrierless, 
highly  exothermic  cases  such  as  this),  the  rate  constant  for 
dissociation  scales  much  more  sharply  with  respect  to  energy 
than  does  that  for  detachment,  with  a  larger  asymptotic  value 
as  well  (see  Fig.  5  in  Ref.  32).  This  is  due  to  the  larger  number 
of  contributing  states  in  the  case  of  dissociation,  since  detach¬ 
ment  is  usually  s-wave,  reducing  the  number  of  potential  rota¬ 
tional  states  that  can  conserve  orbital  angular  momentum.  To 
test  the  applicability  an  approach  similar  to  that  employed  in 
Ref.  32  has  been  attempted  for  the  current  system,  calculating 
the  specific  rate  curve  for  dissociation  of  J  =  0  NCF”*  to  O” 
+  NO  using  Orbiting  Transition  State-Phase  Space  Theory 
and  that  of  electron  detachment  via  microscopic  reversibility 
by  fitting  to  experimental  data  of  thermal  associative  electron 
attachment  to  NO2.  At  the  energy  of  the  reactants,  the  ratio 
of  the  resulting  calculated  number  of  states  for  dissociation  to 
that  of  detachment  is  20:1,  suggesting  a  branching  of  >95%, 
consistent  with  our  experimental  results.  However,  the  calcu¬ 
lated  absolute  rate  constants  approach  1014  s1 ,  a  timescale 
competitive  with  that  of  intramolecular  vibrational  energy 


redistribution,  suggesting  that  a  statistical  approach  may  not 
be  appropriate.  As  such,  a  quantitative  description  of  the  reac¬ 
tion  appears  to  require  a  complete  PES  description,  which  is 
not  currently  available,  and  we  are  limited  to  the  above  quali¬ 
tative  discussion  of  the  branching. 

CONCLUSIONS 

In  conclusion,  we  have  reported  the  rate  constants  for  the 
reactions  of  O'  and  C>2_  with  N  and  O  atoms  from  173  to 
500  K  for  the  reactions  with  O  and  173-400  K  for  the  reac¬ 
tions  with  O2  for  the  first  time.  The  room  temperature  values 
agreed  reasonably  well  with  previous  measurements,  with  the 
exception  of  the  reaction  of  O 2 _  +  O.  The  discrepancy  in 
this  case  is  likely  due  to  the  absence  of  significant  signal  drift 
arising  from  surface  oxidation  instilling  insulating  layers,  an 
effect  minimized  due  to  the  oil-free  nature  of  the  current  sys¬ 
tem.  Additionally,  in  contrast  to  previous  reports,  we  have 
found  the  C>2_  +  N  reacts  to  form  predominantly  (>85%) 
NO  +  O-,  at  all  temperatures  where  a  branching  determi¬ 
nation  was  available.  We  attribute  this  to  specific  rate  curves 
for  detachment  rising  more  slowly  with  energy  than  for  dis¬ 
sociation.  Importantly,  all  of  the  studied  reactions  were  found 
to  display  negative  temperature  dependences,  increasing  their 
importance  in  mesospheric  chemistry,  where  the  lower  tem¬ 
peratures  ensure  that  these  reactions  will  be  occurring  faster 
than  the  previously  reported  room  temperature  rates. 
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